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acceleration  of  gravity,  ft/sec 
vehicle  vertical  altitude,  ft 
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earth  surface  velocity  at  gun  muzzle,  fps 

angle  between  vectors  V  and  b^,,  radians 
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non-dimensional  altitude  increase 
parameters  introduced  in  solutions  of  integrals 
range  increment  angle,  radians 


.  0  INTRODUCTION 


The  precise  calculation  of  the  motion  of  a  gun-launched  vehicle 
through  and  out  of  the  atmosphere  depends  on  a  large  number  of  parameters. 
These  include  the  environmental  parameters  of  atmospheric  properties  and 
gravity,  which  vary  with  altitude,  and  effects  of  the  earth's  rotation, 
which  depend  also  on  the  launch  direction.  Then  there  are  the  vehicle 
parameters  of  mass,  external  size  and  shape,  notion  of  control  surfaces, 
specific  impulse,  mass  ratio,  and  ignition  and  burning  time  of  any  rocket 
motors,  and  auxiliary  or  directional  thrust  systems.  Finally  there  are 
the  gun  launch  parameters  of  initial  velocity  and  direction. 

Inclusion  of  the  accurate  variation  of  these  parameters  along  the 
trajectory  leads  to  equations  of  motion  which  cannot  be  solved  exactly. 
Consequently,  various  numerical  computer  programs  have  been  devised,  and 
these  satisfactorily  determine  trajectories  for  given  values  of  the  para” 
meters.  In  general  they  also  can  rapidly  display  the  results  of  arbitrary 
variation  of  the  parameters,  so  that  optimum  trajectories  among  those 
calculated  can  be  selected. 

This  approach  to  trajectory  calculation,  while  probably  essential 
for  accurate  final  calculations  for  a  particular  vehicle  missior,  suffers 
from  two  serious  defects  when  appLied  to  preliminary  trajectory  calculations, 
in  which  possible  missions  for  existing  or  proposed  vehicles  are  icing  con¬ 
sidered.  First,  a  computer  may  not  oe  as  readily  available  as  the  require- 
merits  for  a  set  of  rapid  preliminary  calculations  would  indicate.  Second, 


and  more  serious,  as  in  all  numerical  programs  the  output  of  numbers  has  to 
be  interpreted  and  analysed  for  trends  in  the  light  of  the  parameter  input. 
There  are  no  analytical  forms  linking  input  and  output,  from  which  trends 
could  be  discerned  and  predictions  made,  and  as  a  result  it  is  difficult 
to  reach  elearcut  decisions  about  new  vehicle  missions. 

It  would  therefore  be  useful  if  simplified  forms  of  the  governing 
equations  could  be  developed  which  would  retain  the  essential  features  of 
the  exact,  equations  and  produce  vehicle  trajectories  and  other  performance 
characteristics  accurate  to  within  5  or  even  10  percent,  while  providing 
straightforward  analytical  or  graphical  links  between  the  various  parameters 
and  the  resulting  performance.  That  is  the  main  purpose  of  the  present 
report.  An  additional  purpose  is  to  provide,  for  reference,  derivations 
of  some  of  the  important  equations  governing  the  motion  of  gun-launched 


vehicles. 


2 . 0  FUNDAMENTAL  ASSUMPTIONS 


It  will  be  assumed  that  the  vehicles  considered  are  aerodynamically 
symmetric  with  respect  to  their  trajectory  at  all  times,  so  that  they 
experience  no  lift,  merely  drag,  in  their  motion  through  the  atmosphere. 

Any  rocket  thrust  will  be  assumed  to  be  constant  in  magnitude,  and  directed 
back  along  the  trajectory,  so  that  the  vehicles  are  turned  only  by  gravity. 
During  motion  of  the  vehicle  either  through  the  atmosphere  or  under  the 
action  of  rocket  thrust,  or  both,  it  will  be  assumed  that  the  gravitational 
force  is  constant  in  magnitude  and  direction.  The  above  assumptions  will 
be  accurate  to  within  3  percent  for  almost  all  cases. 

2.1  BASIC  EQUATION  OF  MOTION  ALONG  TRAJECTORY 

In  the  light  of  the  above  assumptions,  Figure  1  serves  to  define, 
the  equation  of  vehicle  motion  along  its  trajectory.  Newton's  Second  Law 
of  Motion  requires  'Resultant  External  Force  along  Trajectory  =  Time  Rate 
of  Change  of  System  Momentum  along  Trajectory'.  The  external  forces  are 
the  tangential  component  of  gravitational  force  mg  sin  0  and  the  aerodynamic 
drag  D.  The  rocket  motor  thrust  F  is  an  internal  force  accounted  for  by 
the  system  momentum  change.  Consider  the  vehicle  at  time  t  when  its  mass 
is  m  and  its  velocity  V.  After  an  infinitesimal  time  increment  At:  the 
rocket  exhaust  of  constant  velocity  Vj  has  expelled  mass  -Am  (the  minus 
sign  preserves  the  algebraic  sign  convention)  and  the  velocity  of  the 
vehicle  has  increased  to  V  +  Av.  Thus  Newton's  Second  Law  becomes: 


-  5 


D  -  mg  sin  9  =  Lim  — 
&t-»0 


■[(m  4-  Am)  (V  4-  AV)  -  Am  (V  -  V.)j-  -{mv| 


At 


dv  4.  u 

=  hit*  +  V .  -r- 
dt  j  dt 


or 


F  -  D  -  mg  sin  9  =  m 


d V 
dt 


(2.1) 


dm 

where  F  =  -  V.  —  =  V.b,  and  b  is  the  constant  motor  burning  rate. 

J  ut  J 

If  Eq.  (2.1)  is  multiplied  through  by  ~  and  then  integrated  from 

an  initial  state  denoted  (  )  ,  on  the  trajectory  to  a  final  state  denoted 

n-1 

(  )  the  result,  using  the  fact  that  the  rate  of  increase  of  altitude, 

n 


dh  _  „  .  n 

—  =  V  sin  9 
dt 


(2.2) 
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V  *»  V  +  V .  In 
n  n-1  j  m 


n- 1 


J 


mV  sin  9 


dh  -  g 


1  i 

n-1 


_,h 

C  p- 


n-1 


dh 
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(2.3) 


or 


&v  =&v  -  Av  -Av 

n  F  D  g 
n  n 


In  Eq.  (2.3),  the  integrals  for  the  velocity  decrements  from  aerodynamic  drag 

and  gravity  cannot  be  evaluated  by  quadratures  without  the  introduction  of 

further  assumptions.  The  integral  for  is  considered  first. 

n 

2.2  AERODYNAMIC  DRAG  OF  GUN-LAUNCHED  SYSTEMS 

Aerodynamic  drag  is  conventionally  expressed  in  terms  of  a  drag 
coefficient  C^  through  the  defining  equation 

D  =  CD  q  A  (2. A) 


where  A  =  reference  area  =  projected  frontal  area  of  body  for  a  rocket  or 


gun-launched  vehicle. 


q  *  dynamic  pressure. 


In  high  speed  aerodynamics  the  drag  is  caused  by  both  the  compres¬ 
sibility  and  friction  of  the  air,  and  therefore  the  drag  coefficient  depends 
upon  both  Mach  number  M  and  Reynolds  number  NR, 


cD  =  eD  (M,  Nr)  (2.5) 


where  M  =  —  =  Vehicle  Mach  number 
a 

N„-  =  Vehicle  Reynolds  number 

*  V 

a  =  local  atmospheric  speed  of  sound 


d  =  vehicle  base  diameter 


=  local  atmospheric  kinematic  viscosity. 


It  is  convenient  to  express  the  dynamic  pressure  in  terms  of  Mach  number. 


q  =  .*EM2 
2 


(2.6) 


where  ^  =  ratio  of  atmospheric  specific  heats  at  constant  pressure  and 
constant  volume 
=  1.40,  assumed  constant 
p  *  local  atmospheric  pressure. 


The  dependence  of  on  M  and  NR  is  complex,  even  for  simple  vehicle 
shapes,  and  it  has  been  common  to  make  approximations  in  performance  analysis 
In  the  present  ana?.ysis,  the  nature  of  gun-launching  permits  a  very  simple 
and  useful  approximation  to  be  made.  All  gun-launched  vehicles  have  initial 


Mach  numbers  greater  than  3,  and  even  without  rocket  boost,  vehicle  Mach 
numbers  will  rarely  fall  much  below  3  during  motion  through  atmosphere 
dense  enough  to  produce  significant  drag.  Therefore  it  is  not  necessary 
to  consider  the  drag  variation  of  the  vehicles  in  subsonic  and  transonic 
flight  (as  is  required  for  conventional  rocket-launched  systems).  Only 
the  drag  variation  in  supersonic  and  hypersonic  motion  need  be  considered 
and  this  is  much  simpler. 


Over  the  range  of  velocities  through  the  atmosphere  experienced 
by  a  vehicle  on  a  typical  gun-launched  mission,  the  variation  of  with 
will  be  quite  small,  and  it  is  assumed  henceforth  that 

CD«CD(M)  (2.7) 

only  for  a  particular  vehicle.  This  variation  can  be  approximated  quite 
accurately  in  the  supersonic  and  hypersonic  range  by  the  relation 


where  K  is  a  constant  chosen  to  provide  the  best  fit  of  Eq.  (2.8)  to  the 
available  data  for  the  vehicle  over  the  Mach  number  range  of  interest.  In 
Figure  2  such  a  fit  is  shown  for  the  Martlet  2A  glide  vehicle,  the  data 
being  obtained  from  Reference  1.  In  the  figure  the  maximum  deviation  of 
the  approximate  from  the  actual  curve  is  7.7  percent,  and  since  it  is  clear 
that  positive  and  negative  deviations  will  tend  to  cancel  in  the  drag 
integral  of  Eq.  (2.3),  and  that  in  any  case  the  actual  values  of  CD  may  not 
be  accurately  known  at  the  time  of  the  preliminary  calculations  for  which 
t  its  analysis  is  intended,  the  approximation  of  by  Eq.  (2.8)  is  seen  to 
be  satisfactory. 


Therefore,  Eq.  (2.4)  becomes 


and  the  drag  integral  can  be  written 


(JL) 


dh 

m  sin  Q 


(2.9) 


(2.10) 


In  the  integrand  in  Eq.  (2.10),  (p/a)  is  a  property  of  the  local  atmosphere, 
and  is  a  function  of  h  only  for  a  given  model  of  the  atmosphere.  As  the 
vehicle  climbs  through  the  atmosphere  after  gun  launch,  sin  0  is  a  slowly 
decreasing  function  of  h  which  of  course  is  unknown  until  the  trajectory 
is  determined.  However,  the  initial  value  is  known,  and  because  of  the  slow 
decrease  with  increasing  h,  a  sufficiently  accurate  average  value  sin  6^ 
can  be  determined  from  the  equivalent  zero-drag  trajectory  between  the  same 


altitude  limits.  The  drag  integral  can  then  be  written 


(2.11) 


For  vehicles  without  rocket  motors,  or  for  glide  sections  of  any 
gun-launched  vehicle  trajectory,  m  is  constant  and  Eq.  (2.11)  can  be  evalu¬ 
ated  directly.  During  motion  through  the  atmosphere  under  rocket  thrust, 
m  decreases  at  a  constant  time  rate,  and  again  the  equivalent  zero-drag 
trajectory  can  be  used  to  approximate  m  as  an  integrable  function  of  h 
between  the  same  altitude  limits. 


The  significant  observation  to  be  made  from  Eq.  (2.11)  is  that, 
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within  the  accuracy  of  Eq.  (2.8),  the  decrement  in  velocity  experienced  by  a 
supersonic  or  hypersonic  gun-launcheo  vehicle  due  to  atmospheric  drag  is  not 
increased  by  increasing  the  launch  velocity  of  the  vehicle,  as  one  might 
expect  intuitively. 

On  the  contrary,  for  a  vehicle  without  rocket  motor,  or  during  a  glide 
section  of  any  vehicle  trajectory,  the  drag  velocity  decrement  is  nearly 
independent  of  the  vehicle  velocity,  which  enters  only  through  its  rather 
small  effect  on  the  value  of  sin  9n  in  Eq.  (2.11),  and  here  the  effect  of 
increased  launch  velocity  is  to  increase  sin  @n  and  thus  reduce  the  drag 
decrement. 

For  a  rocket-powered  section  of  vehicle  trajectory  with  a  motor  of 
given  burning  rate,  the  effect  of  increased  launch  velocity  is  to  reduce 
the  decrease  of  vehicle  mass  m  in  passing  through  a  given  altitude  increment, 
and  thus  again  reduce  the  drag  velocity  decrement. 

In  general,  then,  Eq.  (2.11)  shows  that  for  any  gun-launched  vehicle, 
the  higher  the  launch  velocity,  the  lower  will  be  the  resulting  decrement  in 
velocity  due  to  aerodynamic  drag. 

2.3  THE  MODEL  ATMOSPHERE 

Before  Eq.  (2.11)  can  be  integrated,  the  dependence  of  (p/a)  on  h  for 
the  atmosphere  under  consideration  must  be  put  in  suitable  analytic  form. 

Here  the  model  used  assumes  a  linear  variation  of  temperature  T  with  altitude, 
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T  =  Vi  -  |3ai  -  Vi)  <2-12> 

where 

£  -  lapse  rate  =  constant. 

The  air  is  assumed  to  obey  the  equation  of  state  of  a  perfect  gas 

p  =  ?  RT  (2.13) 

where  o  -  air  density 

R  =  specific  gas  constant 
=  1716  ft2/ (sec2.  °Rankine)~1. 

The  additional  relation  needed  is  the  vertical  equilibrium  equation  for 
the  static  atmosphere 

a?  -  -  e  *  <2-14> 

If  P  and  T  are  eliminated  from  Eqs,  (2.12),  (2.13),  and  (2.14)  the 
result  is 
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The  speed  of  sound  in  a  perfect  gas  is  given  by 

a  *  /0f  RT 

(2.16) 

and  so,  using  Eq.  (2.12), 

s  ■  vi>} 

(2.17) 

Therefore,  the  required  function  for  (p/a)  is,  with  a  little  rearrangement 


T  -  0(h  -  h  ) 

n-1  v  n“-l 

'  a  ,  / 
n-1 

-  ^h“l 

(2.18) 


For  a  particular  atmosphere,  (p  is  chcsen  to  fit  the  actual  temperature 
variation  within  given  altitude  limits.  In  this  report,  the  reference 
atmosphere  is  the  Cape  Kennedy  Standard  Atmosphere  from  Reference  2,  and 
a  very  good  fit  is  obtained  using  3  values  of  ^  ,  as  follows: 

0<h  <  50.8  kft  -  3.54oR/kft 

50.8  <h  <162.4  kft  ^  =  -1.235°R/kft 

162.4  <h  <  270.6  kft  ^  =  1.574°R/kft 

Figures  3  and  4  show  the  close  agreement  between  the  pressure  and  sound 
speed  from  Reference  2  and  the  values  given  by  Eqs.  (2.15)  and  (2.17), 
using  the  above  values  of  ^  .  For  this  atmosphere,  at  sea  level  (h  =  0) : 
pQ  =  2125  psf,  aQ  =  1138  fps,  g  =  32.15  ft/sec^ 
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3.0  GLIDE  TRAJECTORIES 


When  a  gun-launched  vehicle  is  not  moving  under  rocket  motor  thrust, 


mass  m  is  constant,  and  Eq.  (2.11)  becomes 


Av  ._ES*_ 

2m  rnrsr 


i 


p  \dh  =  _ 

a '  2m  s in  9. 


(3.1) 


The  integral  A  In  is  evaluated  using  Eq.  (2.18)  and  the  solution  of  Eq.  (2.3) 
can  then  be  considered. 


3.1  EVALUATION  OF  DRAG  INTEGRAL 


Using  the  parameters  for  the  model  of  the  Cape  Kennedy  Standard  Atmosphere, 
this  function  has  been  calculated  for  increments  of  altitude  from  sea  level 
until  incremental  values  of  the  function  become  negligible.  The  results  are 
given  in  Table  1,  and  the  total  function  1^  from  that  Table  is  plotted  in 
Figure  5.  It  can  be  seen  that  the  drag  decrement  in  vehicle  velocity  becomes 
negligible  above  160  kft. 


C/1PE  KEm£DV  fimosPHEPE 
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3.2  ZERO-DRAG  GLIDE  TRAJECTORIES 


The  equations  of  motion  for  a  vehicle  moving  without  rocket  thrust 
or  atmospheric  drag  under  constant  g  are  very  simple.  They  are  useful 
in  the  present  analysis  both  in  providing  a  basis  for  approximating  certair. 
functions  in  Eq.  (2.3),  and  in  themselves  as  giving  sufficiently  accurate 
results  for  the  upper  parts  of  vehicle  trajectories  high  enough  to  neglect 
aerodynamic  drag  but  not  so  high  as  to  require  the  inclusion  of  g  variations. 


The  first  relation  needed  is  obtained  directly  from  the  conservation 
of  mechanical  potential  and  kinetic  energy: 


h  -  h 
n  n-1 


V  2  .  u  V  2 
n-1  n 


(3.3) 


The  other  relations  needed  are  obtained  by  manipulating  the  equations  of 
vehicle  motion  in  the  X-  and  h-  directions  (see  Figure  1),  and  are  derived 


in  Appendix  A: 


6K  ■  tan  VlAXn  “  '%  se°2  °-4) 


sin  0  =A - 

n  V  1  -  Ak 


* 2  9n-l  r~  Cos20n-1 

l  - l-A\ 


(3.5) 


where 


Cos  0 


Cos  0n  = 


-V1  "  A?un 


A&  «  (h  -  h  ),  (X  -  x  ) 

n  yZ  n  n-1  n  y-  n  n-1 


(3.6) 


(3.7) 
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3,3  APPROXIMATIONS  FOR  GLIDE  TRAJECTORIES  WITH  DRAG 

For  a  vehicle  moving  without  rocket  thrust,  Eq.  (2.3)  can  be 
written,  using  Eq.  (3.1): 


Avn=-Av  -Av 

n  6n 


-  —  Ai  - 


2m  sin  y 


n 


dh 
*  V 


(3.8) 


n 


n-1 


If  there  were  no  aerodynamic  drag,  Eq.  (3.3)  could  be  used,  and  written 
in  the  form  showing  the  velocity  decrement  caused  by  gravity: 

2S  <hn  -  W 


V  -  V  , 
n  n-1 


Vn  +  V  . 
n  n-1 


or 


where 


Av  =  -  =  Av 

n  -  gn 

n 


V  =  %  (V  ,  +  V  )  -  V  +  %  Av 
n  n"l  n  n-1  n 


(3.9) 


(3.10) 


It  is  now  assumed  chat AVgn  in  Eq.  (3.8)  can  be  approximated  by  an  expres¬ 
sion  of  the  form  of  Eq.  (3.9),  with  Vn  the  arithmetic  mean  of  the  actual 


V  .  and  V  : 
n-1  n 


.  ,  g  Ah 

Av  =  -  Avn  - 1 — a 

n  11  — 


n  V. 


n 


or,  multiplying  through  by  Vn  and  using  Eq.  (3.10): 


%  Avn  +  (Vn-1  +  h  Avn  +  (vn-1AvDn  +  gAV  =  o. 
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Solving  for  A  V^}  and  dividing  through  by 

AY 


&v  /  Av  Avd 

v  -  *  -  a  +  4  ^  >  +  V<i**r*>*-*r*  • 

n-1  n-1  n-1  n-1 


n 


or 


where 


Air  -  -  (1  +  %Air  )  +  /  (1  -  %A>  o  >2  -  A&n 

n  V  n 


Ayn 

Air„  =  —  ■  Air 


AV. 


n-1 


D  V  . 
n  n-1 


(3.11) 


(3.12) 


Before  Eq.  (3.11)  car  be  used  to  calculate  velocity  decrements  along  a 
glide  trajectory,  a  method  of  determining  sin  9^  in  the  expression  for 

AVp  must  be  given.  It  can  be  shown  that 

n  u 

f  h 

1  1  n 

1  1  dh 


sin  9 


n 


is  a  good  approximation  to 


r 

I 


n-1 


n-1 


dh 


sin  9 


where 


sin  9  =  K  (sin  9  .  +  sin  9  ) 

n  z  v  n-1  n 


(3.13) 


and  sin  9^  is  given  by  Eq.  (3.5). 

For  example,  with  sin  9  =  0.50  and  sin  9^  =  0.25,  the  approximate  formula 


is  only  1%  higher  than  the  exact  integral. 


Therefore  sin  9  as  defined  by  Eq.  (3.13)  is  used  in  Eq.  v3.1) 
n 

In  the  calculation  procedure  next  to  be  described,  the  fundamental  inde¬ 
pendent  variable  is  the  altitude  increment  A  h^.  Calculation  of  the 
vehicle  trajectory  requires  corresponding  values  of  the  horizontal 
range  increment  A  X^,  shown  in  Figure  6.  Since 

Ax  =  cot  9*1  Ah 

n  n  n 

Ay  — 

and  8  will  not  b>,  appreciably  different  from  9  as  given  by  Eq.  (3.13) 
n  n 

it  is  assumed  that  a  satisfactory  approximation  is 


Ax  =  cot  9  Ah 
n  n  n 


In  the  same  way  it  is  assumed  that  the  time  increment  A tn  is  given  by 


At  = 


Ah 


n  V  sin  9 
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3.4  CALCULATION  PROCEDURE  FOR  CONSTANT-g  GLIDE  TRAJECTORIES 

For  a  given  vehicle  launched  under  given  conditions  in  a  given 
model  atmosphere,  the  following  parameters  are  known: 

g,  m,  A,  K,  ,  Vq,  6q 

The  factor  in  &V_  can  therefore  be  calculated.  The  drag  integral  Al 
2m  D  n. 

n 

can  be  determined  for  each  altitude  increment  Ah^  from  Figure  5  or  Table 
1.  The  altitude  to  which  the  aerodynamic  drag  is  significant  is  broken 
down  into  a  few  increments,  the  number  increasing  with  the  accuracy  desired. 
Thus  four  increments  should  be  sufficient  for  most  preliminary  calculations, 
while  one  may  be  enough  for  a  rough  calculation. 

The  procedure  is  then  to  use  the  equations  of  3.2  and  3.3  to 
determine  the  velocity  decrement,  horizontal  range  increment,  time  increment, 
and  final  elevation  angle  for  each  altitude  increment.  Above  the  final 
atmospheric  altitude  increment,  the  equations  of  jf  3.2  and  Appendix  A 
are  used  to  calculate  the  remainder  of  the  trajectory  to  its  apogee. 
Alternatively,  Eq.  (3.4)  for  the  parabolic  zero-drag  trajectory  can  be  put 
in  the  universal  form: 

>1  =  2  ^  -  f  2  (3.16) 

where  ^  and  ^  are  defined  in  Appendix  A.  Eq.  (3.16)  is  plotted  in 
Figure  7  and  all  zero-drag,  constant-g  glide  trajectories  can  be  determined 


from  it. 


Figure  7  Universal  Zero-Drag  Constant-g  Trajectory 
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In  greater  detail,  the  procedure  for  the  atmospheric  part  of  the 

trajectory,  given  h  , ,  Ah  ,  V  .,9  , ,  is  to  calculate  A$\,  from 
J  n-1  n  n-1  n-1  n 

Eq.  (3.7)  and  use  this  to  calculate  sin  9  from  Eq.  (3.5)  or  cos9  from 

n  ^  n 

Eq.  (3.6).  (The  latter  is  mora  convenient  for  values  of  9^  near  90°. ) 

This  determines  the  new  elevation  angle  9  .  Next  sin  9  is  found  from 

n  n 

Eq.  (3.13),  giving  9^.  With  Alfi  determined  for  the  given  A  b  from 

Figure  5  or  Table  1,  A  is  calculated  from  Eq.  (3.8),  and  A"lTD  from 

n  n 

Eq.  (3.12).  Now  A"lTn  can  be  calculated  from  Eq.  (3.11),  and  this  gives  the 

velocity  change  Av^  from  Eq.  (3.12).  The  range  increment  A  X^  is  calculated 

from  Eq.  (3.14)  and  this  gives  the  new  corresponding  to  h^.  The  time 

increment  A t^  is  calculated  from  Eq.  (3.15),  giving  the  new  t  .  The  new 

velocity  V  is  given  by  adding  Av  to  V  , ,  The  procedure  is  then 
n  n  n-1  r 

repeated  for  the  next  altitude  increment. 

Figure  8  gives  the  application  of  the  procedure  to  the  trajectory 
of  the  Martlet  2A,  IOWA  shot  from  the  Barbados  16"  gun  on  March  1965. 
Results  of  calculations  with  both  four  and  one  atmospheric  altitude  increments 
are  compared  with  radar-determined  trajectories  of  the  actual  shot  from 
Reference  1,  and  with  the  standard  MRP  trajectory  from  Reference  3.  Carves 
are  given  for  vehicle  weights  of  both  170  lbs  and  180  lbs,  since  actual 
vehicle  weighed  170  lbs,  while  the  HARP  trajectory  was  calculated  for  180  lbs. 
Figures  9  and  10  show  the  variation  of  vehicle  velocity  and  elapsed  time  with 
altitude,  up  to  the  apogee,  for  the  180  lb  vehicle  as  calculated  by  the  present 
method  with  both  four  and  one  atmospheric  steps,  and  by  the  standard  HARP 
computer  program.  Calculations  for  Figures  8,  9  and  10  are  tabulated  in 
Appendix  B.  In  Figure  11,  the  180  lb  trajectories  a' e  compared  again  for 
greater  clarity. 
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Figure  II  Martlet  2A  Trajectories 
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A  comparison  that  makes  greater  demands  on  the  method  is  shown  in 
Figure  12.  Here  a  projectile  is  analyzed  from  launch  to  impact,  with  the 
entire  trajectory  low  enough  that  measurable  aerodynamic  drag  is  experienced 
throughout.  The  present  method  with  four  atmospheric  altitude  increments 
for  both  climb  to  and  descent  from  80,000  ft  and  altitudes  above  80,000  ft 
assumed  non-atmospheric,  is  compared  with  the  standard  HARP  computer  program, 
and  elapsed  time,  velocity,  and  altitude  are  plotted  against  range.  Cal¬ 
culations  are  given  in  Appendix  C. 

3.5  DISCUSSION 


The  calculation  procedure  is  very  simple  and,  even  for  a  choice 
of  several  altitude  increments  in  the  atmosphere,  it  can  be  carried  out 
quickly  using  just  a  slide  rule  or  a  slide  rule  and  trigonometric  tables. 
Despite  this  simplicity,  the  curves  of  Figures  8,  9,  10,  11  and  12  show 
that  it  produces  accurate  results,  agreeing  closely  both  with  actual 
measurements  and  with  more  exact  numerical  procedures. 

Figure  8  shows  that  a  four-step  atmospheric  calculation  predicts 
the  observed  trajectory  of  Martlet  2A  IOWA  very  closely,  with  the  apogee 
only  1.8%  higher  than  that  observed. 

Figures  9,  10  and  11  show  that  the  four-step  atmospheric  cal¬ 
culation  for  a  180  lb  Martlet  2A  gives  results  nearly  identical  with  those 
of  the  HARP  computer  program.  Moreover,  even  the  one-step  atmospheric  cal¬ 
culation  gives  good  accuracy  in  this  example,  the  velocity  values  and 
elapsed  times  lying  quite  close  to  the  HARP  curves,  and  the  trajectory 
having  an  apogee  only  3.0%  high. 
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Figure  12  shows  that  even  a  low-altitude  projectile  trajectory  is 
calculated  from  launch  to  impact  with  good  accuracy  by  the  method.  Errors 
in  calculated  quantities  at  points  along  the  trajectory  are  generally  less 
than  5%,  and  in  some  of  the  more  important  quantities  the  errors  are: 
range  (+1.9%),  apogee  (+4.2%),  range  at  apogee  (+1.1%),  velocity  at  apogee 
(-3.5%)?  time  to  impact  (+4.1%),  time  to  apogee  (+3.3%),  velocity  at  impact 
(-9.0%,).  The  larger  error  in  impact  velocity  is  caused  by  matching  the 
formula  for  drag  coefficient,  =  —  ,  to  the  curve  used  in  the  computer 
program  at  M  =  4.4.  This  gives  excellent  agreement  at  the  higher  Mach 
numbers,  but  produces  higher  drag  at  the  lower  Mach  numbers  preceding  impact. 
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4.0  ROCKET  POWERED  TRAJECTORIES 

The  altitude  change  during  the  burning  of  a  rocket  motor  stage  of 

a  vehicle  is  in  general  small  enough  to  permit  neglect  of  the  change  in  g 

during  motor  burning.  For  altitudes  below  about  160  kft,  however,  effects 

of  aerodynamic  drag  cannot  be  neglected,  and  the  vehicle  velocity  increment 

V  resulting  from  rocket  thrust  is  then  given  by  Eq.  (2.3),  with  the  drag 

velocity  decrement  V  given  by  Eq.  (2.11)  and  the  gravity  velocity 

n 

decrement  A.V  given  by  Eqs.  (3.9)  and  (3.10).  In  this  section  the 

gn 

evaluation  of  &V  from  Eq.  (2.3)  and  the  determination  of  other  rocket- 
n 

powered  trajectory  parameters  are  described. 

4.1  EVALUATION  OF  ELEVATION  ANGLE 

Before  A.VD  can  be  evaluated  from  Eq.  (2.11),  a  suitable 
n  _ 

determination  of  the  mean  value  sin  0^  must  be  made  for  a  rocket-powered 
trajectory,  and  the  dependence  of  mass  m  on  altitude  h  during  thrusting 
must  be  expressed  so  that  the  integral  can  be  evaluated. 

As  in  the  glide  trajectories  of  the  previous  section,  the  deter¬ 
mination  of  sin  9^  is  approached  by  neglecting  the  direct  effect  of  drag 
on  elevation  angle.  It  is  shown  in  Appendix  D  that  in  zero-drag,  ecnstant-g 
conditions,  the  differencial  equation  governing  the  variation  of  elevation 
angle  9  with  time  t  under  the  action  of  thrust  F  is 

2 

— 2  +  (2  tan  9  -  cK.  Sec  9)(^|-)“  =  0  (4.1) 


where 


o(  F 

w~ 
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For  the  constant  thrust  rocket  motors  under  consideration  here,  oC  varies 
with  time  and  Eq.  (4.1)  is  a  fairly  complicated  nonlinear  differential 
equation.  Two  special  cases  of  the  equation  can,  however,  be  treated 
easily.  If  cK.  =  0,  of  course,  Eq.  (4.1)  governs  zero-drag  glide  trajectories 
and  the  solution  of  Eqs,  (3.5)  and  (3.6)  is  obtained.  If  oC  =  constant,  ic 
is  shown  in  Appendix  D  that  Eq.  (4.1)  is  easily  reduced  to  the  integral  form: 


where 


cos  ^~2  9  d  9 
(1  +  sin  Q)0^ 


oC-1  .  sy 

cos  0  .  A  L 

_ _ n-1  n 

(1  +  sin  9  .)** 

n-1 


g  AT 


n 


V 

n- 


1 


(4.2) 


This  is  integrable  by  quadratures  for  oC  =  0,  1,  2,  3,  and  it  is  readily 
evaluated  graphically  or  numerically  for  other  values  of  oC  .  Eq.  (4.2) 
shows  that  ©n  =  9n  (9^,  oC  ,  ATn>*  Now,  A^n>  defined  by  Eq.  (3.7)  is 
also  given  by  the  functional  form  =  A*fl  (9  ,  cC  ,  A'i' ) ,  so 


that  0^  can  be  expressed 


e„  -  Sn 


(4.3) 


Eq.  (4.3)  must  reduce  to  the  value  given  by  Eqs.  (3.5)  or  (3.6)  for  cC  =  0, 
and  must  reduce  to  0^  =  0  for  oC~>0O.  It  is  therefore  plausible  to 
try  as  an  approximation  suitable  for  easy  calculation  the  analytic  forms 


sin  9 

n 


sir/  0  -  f  (OC)  Afc,  ' 

_  n-i  n 

1  -  f  (OC)  Airo 

n 


(4.4) 


-yi  -  f  (cc)Arun 


where 


f  (o)  =  1,  f  (00)  =  0. 


An  exponential  form  for  f(cC)  suggests  itself,  and  it  was  determined 
empirically  as  follows.  It  is  shown  in  Appendix  D  that,  for  constant  OC  } 


(CC  +  sin  0  , ) 

n-1 


sin  9  d  l 
c£  +  sin  9 


+  2  (oC2 


6TC 


sin  0 


I  oC  +  sin  0 


Eqs.  (4.2)  and  (4.6)  were  integrated  numerically  to  obtain  exact  solutions 
for  constant  oC  with  which  values  obtained  from  Eq.  (4.4)  with  different 
choices  of  £(cC)  could  be  compared.  The  final  choice  of  f(oC)  was 

....  -0.250CC0*7  „ 

f(oC)  =  e  (4.7) 

In  Fig.  13  this  function  is  compared  with  the  exact  solutions  at  = 

.25  for  =  30°  and  60°  over  the  relevant  range  of  cC.  .  f(cC)  was 

chosen  to  lie  somewhat  below  the  exact  curves  for  two  reasons.  The  first  is 
to  take  some  account  of  the  effect  of  drag,  otherwise  neglected  i.r.  this 
analysis.  The  second  is  to  relate  to  the  actual  constant-thrust  rocket 
motors  under  consideration,  for  which  oC  increases  steadily  from  ignition  h 
burn-out.  In  applying  the  analysis  to  such  motors,  a  mean  value  of  cC  , 


=  h  (oC 


n-l 


+  cC.  ) 

Li 


(4.8) 


is  used  for  simplicity,  but  this  tends  to  make  f(oC  )  too  small,  and  there- 

n 

fore  sin  9^  too  large.  With  sin  9^  given  by  Eq.  (4.4)  and  (4.7)  si"*  9^  is 
given  by  Eq.  (3.13) . 
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4.2  EVALUATION  OF  DRAG  INTEGRAL 


Since  the  drag  velocity  decrement  is  a  small  fraction  of  the 
velocity  achieved  under  rocket  thrust,  the  mass-altitude  relation  needed 
t<?  determine  it  can  be  based  with  sufficient  accuracy  on  the  equivalent 
zero-drag  trajectory,  with  sin  9  =  sin  9^.  Eq.  (2.1)  can  be  written 


dV  ,,  dm  — : — r- 

m  —  =  -  V .  77  -  mg  srn  9 

dt  i  dt  n 

Jn 


(4.9) 


On  multiplication  by  —  and  integration  this  gives 

V  =  — - ~  =  V  -  V  In  -2—  -  g  sin" 9  (t  -  t  )  (4.10) 

.  r->,  dt  n-1  i  m  ,  n  n~  1. 

sin  9  '  n  n-1 

n 


A  second  integration,  using  the  motor  burning  relation 


m=m  ,  -  b  (t-c  , ) 
n-1  n  n-1 


(4.11) 


gives  _  _ 

(V  .  +  V.  )  sin  9  V.  sin  9 

n-1  j  n  j  n 

,  ,  .  n  ,  v  ,  n  ,  m 

h  =  h  .  i - r -  (m  ,  -  m)  + - r - m  1  n  - 

n-1  b  n-1  b  m  , 

n  n  n-1 


g  sin  9 


—  (Vl  -  B) 


(4.12) 


In  Eq.  (4.12)  it  should  be  noted  that  the  last  term  is  very  small  in  comparison, 
with  the  second  term  on  the  right  side.  Thus,  their  ratio  is 


g  sin  9  (ra  -  m) 

_ n  n-1 _ 

2b  (V  +  V .  ) 
n  n-1  l 

n 


g  sin  6n  (c  -  t  ) 

'  2 (V  +  V  .  ) -  =  0rder  (l0) 

n-1 


using  typical  values  for  the  parameters.  Accordingly,  the  last  term  can  be 


neglected  in  the  expression  for  m  for  use  in  Eq.  (2.11).  The  remaining  terms 
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are  made  dimensionless  by  multiplication  by 


m  ,  V  .  s  m  9 
n-1  j  n 

Jn 


giving: 


AH  * 


n 


where 


m  ,  V .  s  m  8 
n-1  j  n 

n 


<h  '  'Vi>  “(ir~  +  ‘V  +  \  ln  H 

Jn 


(4.13) 


v,  -  ^  ,  n  =  jl. 

n  n-1  n-1 


Eq.  (4.13)  is  plotted  for  relevant  values  of  qj*.  and  in  Fig.  14. 

Jn 

(The  subscript  n  of  course  refers  to  the  evaluation  of  the  function  for  an 

altitude  increment  Ah  =  h  -  h  . .) 

n  n  n-1  • 


Although  a  simple  expression,  Eq.  (4.13)  is  not  directly  useful 
in  evaluating  Eq.  (2.11),  since  it  is  transcendental  and  VV|/  cannot  be 
given  explicitly  as  a  function  of  Ah.  However,  is  unity  for  AH  =  0 

and  approaches  zero  as  AH  becomes  large,  so  approximation  by  a 

function  of  the  form 


•wy,  - .  -£CvV  Ah  <4 

is  plausible.  It  was  found  empirically  by  comparison  with  the  curves  of 
Fig.  14  that  the  function 


f (\jv  )  =  i.io  nr. 

I 


0.60 


(4.15) 


gives  satisfactory  agreement  with  Eq.  (4.13).  Thus,  Eq.  (4.14),  using 


/o 
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Unfortunately,  however,  the  function  for  (p/a)  given  by  Eq.  (2.18) 
is  riot  integrable  by  quadratures  in  combination  with  Eq.  (4.14).  It  is 
therefore  replaced  for  this  integration  by  the  function 

(!)  -(f)  e-C<h-V1>  (4.18) 

n-1 

where  c  is  determined  empirically  to  give  good  agreement  with  the  true  vari- 
at  ion  of  (p/a)  over  a  selected  altitude  range.  Eq.  (4.18)  would  be  exact 
for  an  isothermal  atmosphere,  and  is  a  good  approximation  for  an  actual  atmo» 
sphere  over  limited  ranges  of  altitude.  Thus,  in  Fig.  15,  Eq.  (4.18)  is  seen 
to  give  excellent  agreement  with  the  Cape  Kennedy  Standard  Atmosphere  for 
altitudes  between  40  kft  and  120  kft,  with  c  =  0.04808/kft. 


With  the  com.  (.nation  of  Eq.  (4.14)  and  (4.18),  A<-Tn  can  be  evalu¬ 


ated: 


PEF.Z 

/?•/  * 

C  «  0.04808 /left-  _ 


"  40  “ 


AJ  Jl)  f  L  +  _iao 

n  U'n-1  J  e)  ITj  »• 

hn-l  L 


where 


in' 


d  |  it.  0.40  ,,  — — E  V 

-  L  mn-l  Vl  *“  9„J 


1.10  b 

n 

0.40  - 

u.  .  m  .  V  ,  sin  9 
Jn  n-1  n-1  n 


(h  -  hn-l> 


(4.19) 


(4.20) 


4.3  DETERMINATION  OF  VELOCITY  AND  TRAJECIORY  PARAMETERS 


It  is  now  possible  to  consider  the  evaluation  of  the  net  velocity 

increment  Av  arising  from  the  effects  of  rocket  motor  thrust,  aero- 
n  5 

dynamic  drag,  and  gravity: 


Av  =  A v_  -  Avn  -  Av 

n  F  D  g 

n  n  6n 


(2.3) 


Avn-  Avf  -  AVD  -  v  ■  *f 

n  n  n-1  n 


using  Eo4.  (3.9)  and  (3.10).  This  is  made  non-dimensional  by  division  by 


where 


Air 


Av  =  Av  -  - 2—  -AV* 

n  F  2  + AV  D 

n  n  n 


F  =  -V-  in-rv^D  -  — K*A  _ 

n  Jn  n  2m  ,  V  .  sin  9 

n-1  n-1  n 


(4.21) 


(4.22) 


-  41  - 


Solving  Eq.  (4.21)  for  AlT^ 


Airn  =  -  i  + 


k  (tiVF  -Av-jj  )  +.  /  1  +  %  (AiTf  -A\Td  >{  -  A&n  W.23) 

n  n  \j  *•  n  a  J 


Eq.  (4.23)  reduces  to  Eq.  (3.11)  if  Air  p  =  0  (no  rocket  thrust).  The 

n 

procedure  for  evaluating  Eq.  (4.23)  is,  however,  different  from  the  previous 
In  evaluating  Eq.  (3.11),  Ah  was  selected  and  therefore  At  was 


one 


n 


given,  and  the  time  increment  A t  was  part  of  the  solution.  In  the 

present  case,  in  which  the  burning  of  one  rocket  motor  stage  constitutes  a 

trajectory  increment,  the  time  increment  At^  is  known,  as  is  the  mass 

decrement  Am  ,  so  must  be  found  from  the  mass-altitude  relation, 

n  n 

In  Eq.  (4.12)  there  is  no  need  to  neglect  the  small  final  term 

at  this  stage  of  the  calculations,  since  only  simple  algebra  is  involved. 

2g 


Therefore,  on  multiplication  by 


,  Eq.  (4.12)  becomes: 


n-1 


Af>  =  2V.  £  Ah 

n  l  u  n  r, 

Jn 


(1  ->)  )' 

n  n 


<r. 


where 


S. 


gm  ,  sin  9 
_  n-1 _ n_ 

V  !  b 
n-1  n 


(4. 2  A 


and  Ah  is  given  by  Eq.  (4.13)  or  Fig.  14.  The  calculation  procedure  can 
n 

now  be  described.  It  is  convenient  to  introduce  a  parameter  ^  1 


c,  - 


V  b 
n-1  n 


g  m 


-.40' 


n-1 


defined  by: 
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In  terms  of  ?  , , 

n-1 


V.  b 
J  n 


°C  =  k  (c<  +  ©C  )  =  - -  (  — — 

n  n-1  p/  2g  mR  ^ 


.  _ L.  \  =  Jn  £  n-1  n  ,  _1_,  (4.26) 

m  '  2  U  V 

n 


sin  9 

and  <fn-  -g--E' 

°  n-1 

A  complete  rocket-powered  stage  is  calculated  as  one  trajectory 

increment.  The  known  initial  parameters  are  V  , ,  t  , ,  h  .,0  ,  ,  X  ,  g, 

n-1’  n-1’  n-i’  n-1’  n-1  b! 

K,  ,  A,  m i ,  b  ,  c,  V.  ,  (p/a)  .  At  is  known,  and  this  determines 

n-  i  n  j  n  n-1  n 

t  and  A  m  .  IT.  and  Th,  are  calculated  and  AH  is  found  from 
n  n  JR  wn  n 

fig.  14.  £  and  oC  are  calculated,  and  a  trial  value  of  sin  0 

n- 1  n  p. 

is  estimated,  based  on  the  knowledge  of  sin  6  ,  and  oC  .  Using  this, 

n-1  r>  °  ’ 

^  is  calculated  from  Eq.  (4.27)  and  A'i  from  Eq.  (4.24).  f(5c  ) 
n  n  n 

is  then  found  from  Eq.  (4.7)  and  cos  0  is  calculated  from  Eq.  (4.5).  This 

n 

determines  sin  9  and  leads  to  a  new  value  of  sin  0  .  If  this  differs 
n  n 

from  the  estimated  value,  the  process  is  repeated  from  that  point,  using 
the  new  value  of  sin  0r.  Only  the  one  iteration  should  ever  be  required. 
Next  Ah  is  calculated  from  A&-»  ,  giving  h  .  is  calculated 

from  Eq.  (4.20)  and  used  to  calculate  AjR  from  Eq.  (4.19).  This  gives 
Air  n  from  Eq-  (4.22),  and  AV*_  is  then  calculated,  and  AV  is 

u  r  ri 

n  n 

determined  from  Eq,  (4.23).  This  gives  Av  and  V  .  Aange  increment 

n  n 

AXn  can  now  be  calculated  from  Eq.  (3.14)  and  this  gives  range  X  ,  com¬ 
pleting  the  stage  calculation. 


As  an  example,  the  method  was  applied  to  a  Martlet  4,  a  g„r - 
launched  three-stage  rocket,  from  launch  to  burn  out  of  the  second  stage 
(the  remainder  of  the  trajectory  is  high  enough  to  require  consideration  oi 
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the  g-variation) .  It  was  launched  at  6.000  kfps  at  an  elevation  angle  of 
33.0°.  The  first  stage  motor  was  ignited  at  40.0  kft,  and  the  second  stage 
motor  immediately  after  burn-out  of  the  first  stage.  Three  altitude  increments 
were  used  in  calculating  the  glide  trajectory  before  ignition.  In  Table  2 
the  results  of  the  calculations  are  compared  with  the  output  of  the  standard 
HARP  Computer  Program  for  the  same  example  obtained  from  Ref.  4.  Vehicle 
parameters  and  trajectory  calculations,  and  the  output  of  the  Computer  Program, 
are  given  in  Appendix  E. 


TABLE  2 


Martlet  4  Trajectory  Calculations  for  HARP  Case  1046 


V 

n 

kfps 

X 

n 

kft 

h  t 

n  n 

kft  sec 

9 

n 

o 

Launch 

6.000 

0 

0  0 

33.0 

Stage  ]  Ignition 

HARP 

4.822 

67.2 

40.0  14.61 

28.92 

Present 

4.808 

66.4 

40.0  14.67 

28.8 

Error 

-.29% 

-1.34% 

+.41% 

-.42% 

Stage  2  Ignition 

HARP 

12.155 

174.8 

93.9  29.61 

26.02 

Present 

12.21 

172.8 

97.1  29.67 

27.7 

Error 

+  .45% 

-1.14% 

+3.41%  +.  20% 

+6.46% 

Stage  2  Burn-out 

HARF 

18.621 

311.8 

158.7  39.61 

25.35 

Present 

18.84 

306.8 

166.9  39.67 

27.3 

Error 

+1.18% 

-1.67%  +5.16%  +.15% 

+7.69% 
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4.4  DISCUSSION 

Table  2  and  Appendix  E  show  that  rocket-powered  trajectories 
through  the  atmosphere  can  be  accurately  and  easily  calculated  by  slide 
rule  by  the  present  method.  In  the  example  the  velocity  and  range  at 
burn-out  of  Stage  2  were  given  very  accurately,  the  altitude  and  elevation 
angle  less  accurately,  both  being  too  high.  The  time  was  of  course  given 
accurately,  since  it  was  known  except  for  the  initial  glide  trajectory. 

The  high  values  of  h  and  8  ,  although  acceptable,  lead  to 

n  n 

over-estimates  of  apogee,  or  of  orbital  height  after  Stage  3  firing,  and 
this  suggests  further  work  in  improving  the  method  of  determining  sin  0^, 
so  as  to  decrease  its  value  towards  the  correct  one  while  retaining  a 
simple  dependence  on  the  relevant  parameters. 
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5.0  HIGH  ALTITUDE  AND  ORBITAL  TRAJECTORIES 

Rocket -powered  vehicles  will  in  general  reach  very  high  altitudes 
or  go  into  earth  orbits.  The  upper  parts  of  their  trajectories  must  then 
be  calculated  by  recognizing  that  gravitational  force  is  a  central  force 
varying  inversely  as  the  square  of  the  central  distance.  No  aerodynamic 
drag  need  be  considered  unless  such  problems  as  the  gradual  decay  of  orbits 
are  being  studied,  and  any  high  altitude  use  of  rocket  thrust  can  be  treated 
as  a  constant-g  problem  using  an  appropriate  value  of  g,  because  of  the 
short  duration  of  the  firing. 

The  problem  then  becomes  that  of  determining  zero-drag  glide 
trajectories  under  the  action  of  an  inverse-square  central  force.  The 
basic  equations  governing  this  motion  are  simple  and  lead  to  explicit 
solutions  for  trajectory  parameters,  so  there  is  no  need  to  seek  even 
simpler  approximate  solutions.  In  this  section  the  usual  solutions  are 
merely  derived  to  complete  the  set  of  trajectory  equations.  No  examples 
are  worked,  since  the  solutions  are  not  new. 

5.1  EFFECT  OF  EARTH  ROTATION 

In  previous  sections  of  the  report,  trajectory  analysis  was 
given  in  terms  of  vehicle  velocity  V  relative  to  earth.  This  was  suitable 
over  short  distances  for  which  gravity  could  be  considered  a  constant 
vector.  When  the  true  central  force  nature  of  gravity  is  considered,  the 
trajectory  analysis  must  be  based  on  the  vehicle  absolute  velocity  of 
magnitude  U,  or  in  vector  form: 


u! 

E 


(5  A) 
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where  is  the  vector  velocity  of  the  earth's  surface  at  the  gun  muzzle. 

»  i  ^ 

Since  U_  and  V  are  in  different  directions,  it  is  convenient  to  express 

Ci 

relations  between  components.  These  relations,  in  the  context  of  this 
report,  would  be  applied  at  the  point  on  the  vehicle  trajectory  where  the 
constant-g  analysis  of  previous  sections  is  replaced  by  the  analysis  of 
this  section.  The  relations  can  be  obtained  from  Fig.  16.  First,  using 
the  cosine  law: 


U2  =  V2  +  U  2  +  2VU„  cos  X  (5.2) 

E  b 

But,  from  the  right  triangles  of  the  diagram  it  can  be  seen  that 

V  cos  9  cos  =  V  cos  X  (5.3) 

so  that,  cancelling  V  and  substituting  for  cos^in  Eq.  (5.2): 

2  2  2 

U  =  V  +  tL,  +  2VU  cos  0  cos  V  (5.4) 

EE  *  v 


where  is  the  azimuth  angle  of  the  gun.  Eq.  (5.4)  gives  U.  The 

absolute  path  angle  (5  to  the  local  horizontal  is  given  by 


U  sin  &  =  V  sin  9 


15.5) 


The  absolute  azimuth  angle  ^  is  given  by 

U  cos  &  cos  =  V  cos  9  cos  ^ ^ 


(5.6) 
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5.2  FUNDAMENTAL  EQUATIONS  AND  SOLUTIONS 

Because  only  the  central  force  of  gravity  acts  on  the  vehicle 
its  trajectory  remains  in  a  plane  passing  through  the  earth’s  centre,  and 
Fig.  17  defines  the  variables  of  the  trajectory.  The  equations  of  motion 
of  the  vehicle  for  the  r-  and  0-  directions  are  written,  for  unit  mass: 

2 


.2  Art,  2  rP  2 

, .  .  .  dr  /  d0  s  /  E  s 

r-drrectron  — r  -  r  (— r—)  =  -  gp  (~ — ) 

dt2  dt  E  r 

0-direction  r  d  +  2  ~  =  0 

-  at  dt 

dt 


(5.; 

(5.1 


The  second  terms  on  the  left  side  of  the  two  equations  are  the  centri¬ 
petal  and  coriolis  accelerations,  respectively,  and  g^  is  the  value  of  g 
at  the  earth's  surface,  used  in  the  previous  sections  of  the  report. 

If  Eq.  (5.8)  is  multiplied  by  r,  it  can  be  written: 


or 


where 
d0 


2  d0  ^ 

r  - —  =  C,  constant  =  rU  cos  6 
dt 


C  =  r  .  U  .  cos  &  . 

n-1  n-1  n-1 


dt 


can  be  eliminated  from  Eq.  (5.7)  using  Eq.  (5.9): 


(5.5 

(5.1 


d2r  _  C2  _  2 

,  2  3  SE  ^r  ^ 

dt  r 


dr 


,  =  U  sin  6  =  U  .  Therefore 

dt  r 


.2  dU  dU  „3  g_,  r„ 

d  r  _  __r  _  __r  _  C_  °E  E 

,2  dt  r  dr  3  2 

dt  r  r 


(5.1 


(5.: 


Now 


Figure  17  Trajectory  Variables 


1 


Eq.  (5.12)  can  be  integrated  directly  to  give 

"r2-ur2  -?>  f°2<r-.+7>-2%rE2}  ( 

n-1  n-1  [_  n-1  J 

Since  the  altitude  change  of  the  vehicle  is  of  primary  interest,  it  is 


(5.13] 


convenient  to  introduce 


A rv=  AJl.  =  1  _  Isii 


(5.14: 


In  terras  of  Ahj  ,  Eq.  (5.13)  appears  in  the  form 

"«*  *  Ut2  .  -t,  < 


(2  -  A  to) 


where 


Urf  =  U  cos  S 
V 


At  apogee,  U  =  0,  so  that  Eq.  (5.15)  can  be  solved  for  Abb  : 

r  “r 


A/v  = 

ap 


L  n-1 J  n-1 


(5.15; 


(5.16 


Eqs.  (5.15)  and  (5.9)  determine  vehicle  velocity  as  a  function  of  altitude. 
It  is  necessary  to  relate  altitude  to  elapsed  time.  This  can  be  done  using 


dr  _  r^d  Ah/  _  rn-l^  AK/ 

dt  ~  0  r  ,U  /m2„ 

r  n-1  r  (1  -Ah)  U 


(5.17 


Eq.  (5.17)  can  "be  integrated  directly,  using  Eq.  (5.15),  to  give 
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In  order  to  find  the  range  of  the  vehicle,  it  is  necessary  to  know  the 
change  in  0.  This  can  be  found  using  liq.  (5.9): 


dt  dr 


U 

r 


C_ 

r2 


Referring  to  Fig.  ''S  the  vehicle  range  increment  is  the  great  circle 

distance  between  the  point  1  intercepted  at  time  t  at  the  earth's  surface 

by  radius  r  and  the  position  2  at  time  t  of  the  point  3  on  the  earth's 

surface  intercepted  at  time  t  .  by  radius  r  . : 

n-i  n-i 
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Ax  =  r  A  0'  (5.22) 

e 

A0'  can  be  obtained  by  solving  the  spherical  triangle  123.  In  this, 

A0  is  known  from  Eq.  (5.21),  ^ ^  is  given  by  Eq.  (5.6)  and  great  circle 

arc  23  and  angle  231  can  be  found,  since  the  latitude  of  3  and  the  small 
circle  (eastward)  distance  32  travelled  by  point  3  in  the  time  At, 
given  by  Eq.  (5.18),  are  known. 
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\ 

I 

i 


A- 2 


If  0  is  set  equal  to  0R  and  (tn  -  t^)  is  elimated  using  Eq.  (A-2) 
the  result  can  fce  written: 

tan  9n  =  tan  0n.1  -  k  sec2  0n_1  AX 


or 


AX  -  ~2(tan  °n-l  ^  tan  9n> 


sec2  .0 


(A- 7) 


n-1 


Again,  if  Eq.  (3.4)  is  solved  for  the  result  is: 

2(tan  9^  -  tan  6n_^  “  sec^®n-l 


^n  = 


(A- 8) 


sec^  0, 


n-l 


On  comparing  Eqs.  (A-7)  and  (A-8)  it  is  seen  that 


tan 


9n=  Vl  "  8ec2  0n-l  A^n 


(A- 9) 


With  a  little  rearrangement,  using  trigonometric  identities,  this  can 
be  put  in  the  form: 


/sin2  9  -  Ai,n  ' 

sinen'y — t —  ' 


z 

cos  9 


1 


1  - 


n-l 


1  -  Afvr 


or 


cos  Q  = 


cos  9n_i 


(3.5) 


(3.6) 


If  Eq, (A-4)  is  multiplied  through  by  "I  — 7“  the  result  is 

Vl8in  ®n-l 

(  2s^^n  \  2  f  \  / fcn  \ 

VVn-lsin2Vl)  V  Vn-lsin0n-l  )  '  \  Vl8inVl/ 


(A- 10) 


ven:  g  “  32.15  ft/ sec2,  O  =  1.400,  A  =  0.1364  sq.ft.,  K  =  0.920 


42.3  272.3  .368  14.28  33.2  2.03  36.0  81.4 

.23.2  303.2  .554  21.5  40.4  1.46  54.2  99.6 

,38.4  318.4  .682  26.5  45.4  66.8  112.2 

53.9  333.9  1.000  38.8  57.7  .40  97.9  143.3  <-Apogee 


APPENDIX  C 

Projectile  Trajectory  -  Present  Method 
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APPENDIX  D 

Trajectory  Equations  for  Vehicle  under  Rocket  Thrust  with  Zero  Drag,  Constant  g 

Newton's  2^  Law  for  Directions  Parallel  and  Normal  to  Trajectory: 

Parallel:  F  -  mg  sin  9  =  m  —  D-I 

dt 


Normal: 


mg  cos  9  =  mV 


.  d9 


Elimate  V  from  D-l  and  D-2  by  dividing  D-2  by  —  and  differentiating  with 


respect  to  t: 


m  fr  =  "  It  ('"d9/df~')  =  mg  [Sin  9  +  C0S  °  Td§7dF)2}=  F-  mg  sin  9 


dt  \  dP/dt 


Therefore,  on  simplifying  the  a.;ove,  the  differential  equation  for  9(t) 


is  obtained: 

2  2 
+  (2  tan  9  -  «Csec  9)  (— )  =  0 

.it2  ’ 

Solution  for  oC  =  constant 


(4.1) 


Set  6J  = 


Then  Eq.  (4.1)  can  ..e  written 


+  (2  tan  9  -  oC  sec  9)  d9  =  0 


lis  can  ue  integrated  directly  to  give 

cos'*'"^  i  (1  +  sin  9)5C 


cos*"  9  (1  +  sin  9n-]i)': 


This  is  integrated  in  turn  to  give,  after  some  simplification: 
C  11  cos'*  q  d e  cos'*-"  ^  9n_-| 


(1  +  sin  9) 


(1  +  sin  Qn_1)0< 


(4.2) 


constant 


Solution  for  At  for  oC  = 
n 


Divide  D-l  by  m  and  differentiate  with  respect  to  t: 

j2'’.  =  -  e  cos  9  i®  =  &2_.pos£e 
dt  V 


dt2 


D-5 


using  D-2.  Again,  if  D-l  is  solved  for  g  sin  0: 


dV 


•  O  S  U  V 

g  sin  9  =  got-  -  — 


D-6 


or 


g2(l  -  cos29)  =  g2oC2  -  2gt*C—+(  —  \ 

Hf-  v  fit-  / 


2^-2 


dt  v  dt 


D-7 


2  2 

Equating  g  cos  6  from  D-5  and  D-7,  9  is  eliminated: 


„  d2V  2, 2  n  .  0  ,  dV  /dV\2 

=  -  s  -i)  + 


or 


dt  '  dt  • 

2 


£(v£  i) 


D-8 


D-3  can  be  integrated  directly  to  give: 


V 

{£  - 2  •*) 

e: 

>• 

ii 

dV 

can 

oe  eliminated, 

using  D 

dt 

V 

+  sin  9^ 

“  Vn-1 

or 

V 

fot+sin 

Vx] 

=  Vn-1  1  <*+ 

sin  9  J 

+  g(oC“  -l)(t-tn_1) 


g(°^  -l)(t-tn_1) 


dh 


<e£  +  sin  9) 


dt  sin9 


This  gives  directly 

.  .  „  ^  sin  9  dt  ,2  sin  9  (t-t_  i)dt 

dh  =  V„  -i  (c<+  sin  9„  TZTl — • — +  g(o<.  -1)  _ n~ 

n-i  a-1'  (oC+  sin  9)  («+  sin  9) 


D-ll  is  made  non-dimensional  oy  multiplication  with 


la_ 


n-l 


2  and  then 


integrated: 


AT 


( 

=  2(<oC+  sin  9  , ) 

**  n  n- 1 


+  2(0C  -1) 
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D-9 


D-10 


D-ll 


sin  9  dT* 
(oC+  sin  9) 


sin  9  AT  ciTf 
(oC+  sin  9) 
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